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We investigate a lattice version of QED by numerical simu- 
lations. For the renormalized charge and mass we find results 
which are consistent with the renormalized charge vanishing 
in the continuum limit. A detailed study of the relation be- 
tween bare and renormalized quantities reveals that the Lan- 
dau pole lies in a region of parameter space which is made 
inaccessible by spontaneous chiral symmetry breaking. 



I. INTRODUCTION 

QED is the best tested of all quantum field theories. 
But all its success is in the context of perturbation theory. 
It has long been known that there are potential problems 
in the foundations of the theory due to the existence of 
the so-called Landau pole jlj]. In the leading logarithmic 
calculation one finds 
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where e (e^) is the bare (renormalized) charge, m R is the 
renormalized fermion mass, Nf is the number of flavors 
and A is the ultra-violet cutoff. When one attempts to 
send the cutoff to infinity while keeping e R fixed, one 
finds that e diverges at 



A = A; 



(2) 



the location of the Landau pole. The problem can also be 
seen by looking at the gauge invariant part of the photon 
propagator 



D(k) _ 1 
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which has a ghost pole at k 2 — A| . This would mean that 
the entire theory is only applicable for momenta smaller 
than A^. On the other hand, when one keeps e fixed and 
sends the cut-off to infinity, the renormalized charge goes 
to zero, meaning that the theory is trivial. The situation 
in two-loop perturbation theory is much the same. The 
(renormalized) (5 function 



Pr = ra R 



del 
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remains positive for all e R , and the Landau pole is dis- 
placed to lower values 
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QED is not the only theory with a Landau pole prob- 
lem. Every theory which is not asymptotically free suf- 
fers from this problem. While Al ~ 10 227 GeV if only the 
electron is considered, A^ ~ 10 34 GeV in the Standard 
Model. In the Minimal Supersymmetric Standard Model 
(MSSM) A L ~ 10 20 GeV, and in the MSSM with four 
Higgses, which offers a solution to the strong CP prob- 
lem, the Landau pole moves down to Al — 10 17 GeV [|. 
Thus the Landau pole is by no means academic. 

To find a solution to this problem, one must consider a 
non-perturbative formulation of QED. Thus it is natural 
to investigate the problem on the lattice. On the lattice 
the inverse lattice spacing takes over the role of the ultra- 
violet cut-off, a -1 ~ A. Early calculations have shown 
that the non-compact formulation of the theory using 
staggered fermions undergoes a second order chiral phase 
transition at strong coupling ||,Q: 
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The solid line am — 0, m being the bare mass, e 2 > 
e 2 is a line of first order chiral phase transitions, where 
am R , a 3 (ijnp} ^ 0, even though the bare mass is zero. 
The dashed line am = 0, e 2 > e 2 is a line of second 
order phase transitions on which am Rl a 3 (ipip) — 0. A 
meaningful continuum limit can be taken at the tricritical 
point am — 0, e 2 = e 2 , because here we can take a to 
zero while keeping m R fixed. 

To understand the continuum limit of the theory, we 
need to know the renormalized charge as a function of the 
cut-off in the critical region. We have recently computed 
the chiral condensate on large lattices ||. In this letter 
we compute am R and e R in order to understand the fate 
of the Landau pole in QED. 
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FIG. 1. The renormalized mass against the bare coupling 
on 12 4 and 16 4 lattices. 
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FIG. 2. The chiral condensate against the renormalized 
mass on 12 4 and 16 4 lattices. 



II. FERMION MASS 

We obtain the renormalized mass from the fermion 
propagator as outlined in Ref. ||. We are using stag- 
gered fermions which in the continuum limit correspond 
to Nf = 4 flavors of dynamical Dirac fermions. The re- 
sults are shown in Fig. 1. 

Fitting and extrapolating this data it greatly helps that 
the chiral condensate a — a 3 ^^) is a function of am R 
only ||. In Fig. 2 we plot a as a function of anij;. We 
observe that a is well described by the polynomial 



a = 0.6197am fl - 0.321(am«) 3 
-0.040(ami?) 7 , 



0.169(am fl ) £ 
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where the first coefficient is given by the one- loop result. 
This helps because we have already found an equation of 
state (EOS) that describes the a data [0. Combining the 
EOS with the polynomial (|^) gives the curves shown in 
Fig. 1 and the extrapolation to am — 0. (Here we have 
used fit 1 of Ref. ||. Our results do not change quali- 
tatively if we use any of the other fits described there.) 
For 1/e 2 < 1/e 2 chiral symmetry is broken, and even at 
am = the renormalized mass is non-zero. This means 
there is an excluded region shown in white in Fig. 1 (the 
accessible region being shown in gray). 



III. RENORMALIZED CHARGE 

The renormalized charge is obtained from the residue 
of the photon propagator, e 2 R = Z 3 e 2 and Z 3 = 
linifc^o limy^oo D(k). We can compute D(k) on the lat- 
tice, but not at k = 0. The smallest momentum that we 
can reach is 2it/aL, where L is the lattice size (L = 16 



and 12 in our case). To extrapolate to k = we need to 
make a fit to the photon propagator. The k dependence 
of the photon propagator is given by 
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2 = -H(k,mR,L) , 
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where II is the polarization function. In the infinite vol- 
ume limit we then have 
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We have already seen that the non-perturbative II is ac- 
tually very close to the result of one-loop renormalized 
perturbation theory ||. So it is reasonable to make an 
ansatz which is inspired by renormalization group im- 
proved two- loop perturbation theory. In Ref. j7] it is 
shown that to next-to-leading logarithmic order the po- 
larization function can be written 



n = u - 



v 



e 2 U) 
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where U is the one-loop perturbative result, and V the 
two-loop one. The lattice result for U is known ||. For 
V we make the ansatz 



(10) 

i 2 R limits. For 
(/3 2 /2)lna 2 fc 2 , 
<C 1 we should have V ~ 
(/?2/2) lna 2 m|j,. The one-loop result U has these proper- 



This is motivated by the small k 2 and 
a 2 k 2 1 we should have V : 
a 2 k 2 



a 2 m 2 R <§; 
and for 
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ties. We fit this ansatz to a total of 52 photon propaga- 
tors on 16 4 and 12 4 lattices for various values of am, e 2 
in the range 0.005 < am < 0.16, 0.16 < 1/e 2 < 0.22 close 
to the critical point at 1/e 2 = 0.19040(9) [§. A plot for 
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FIG. 3. The residue of the photon propagator against the 
momentum for 1/e 2 = 0.20, am — 0.005 on the 16 4 lattice. 
The open symbols are the data, the solid symbols are the fit. 

one particular parameter set is shown in Fig. 3. For the 
fit parameters we obtain vo = —0.00207(2) and v\ = 
—0.0328(7), giving x 2 /d.o.f. = 1.7. Two-loop continuum 
perturbation theory would give v\ = flz/fti = 3/167T 2 = 
0.0190. In Fig. 4 we show the resulting f3 function for 
e 2 = e 2 . We compare this with the one-loop result. We 
see that the /3 function is a little smaller than the one- 
loop value and is positive. In particular this means that 
there is no ultra-violet stable zero in the /3 function out 
to 4, = e 2 , the maximal value efj can take because Z3 < 
1 |^||. As arriR — > 00 fermion loops are suppressed and 
e R — > e 2 , so that the (3 function vanishes. But this is of 
course not an interesting zero of the j3 function. 
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FIG. 4. The f3 function against the renormalized charge. 
The solid curve is our result, the dashed curve is the lattice 
one- loop result. 



FIG. 5. A sketch of the mapping from the bare parameter 
plane (top) to the renormalized parameter plane (bottom). 



IV. THE LANDAU POLE 

Having calculated the renormalized mass and charge, 
we are now able to discuss the mapping from the bare 
parameters am, e to the renormalized parameters arriR, 
eR. Qualitatively this is displayed in Fig. 5. One can 
choose any e 2 > 0, am shown in the top part of the 
figure by the gray region. This is then mapped onto 
the corresponding gray region in the bottom part of the 
figure. The line am = from e 2 = to e 2 = e 2 is 



For e 2 > e 2 we 



mapped onto the point am^ = e R — 
have already seen in Fig. 1 that am R > 0, even when 
am = 0, because of chiral symmetry breaking. Thus the 
line am = from e 2 = e 2 to e 2 = 00 is mapped onto 
the border line of the gray regions, and the white area 
is inaccessible for any combination of bare parameters. 
From this figure we evidently have triviality. Removing 
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FIG. 6. The (quantitative) mapping of the bare parameter 
plane (top) to the renormalized parameter plane (bottom). 

the cut-off (a — > 0) is only possible at e 2 R = 0. At any 
finite value of e 2 R there is a minimal possible value for 
|am#|, namely the boundary of the gray region. The 
position of the Landau pole is sketched by the dotted 
line. 

We now turn to the quantitative analysis of the prob- 
lem. In Fig. 6 we plot again the bare and renormalized 
planes, this time using 1/e 2 and 1/efj, respectively, as the 
horizontal variables because this displays the asymptotic 
behavior best. The curves are lines of constant l/e R 
(top part) and lines of constant am (bottom part), re- 
spectively. All lines of constant e R end on the first order 
phase transition line, and only the line e R = goes into 
the critical point. This is another expression of trivial- 
ity of the theory. In the bottom part of the figure the 
gray region is again the allowed region, and the white 
region is inaccessible. The border line is the line am = 0. 



From eq. (g) we find the Landau pole by setting the bare 
charge to infinity. This gives the dotted line. We see 
that it completely lies in the inaccessible region. It runs 
roughly parallel to the border line am = 0. 

We also want to be sure that the photon propagator 
has no ghost pole for any k 2 . Beyond leading logarithmic 
order the extra pole in the photon propagator and the 
divergence in e 2 need not appear at the same place. The 
ghost pole position (in the infinite volume) is given by 
l/k 2 D(k) = 0. It has a solution if 

— < maxLI(/c, m#, oo) — 11(0, trr, oo). (11) 
e R k 

The solution is given by the dashed line in the bottom 
part of Fig. 6. This lies close to the Landau pole, even 
deeper in the inaccessible region. 

V. CONCLUSION 

From Figs. 5, 6 we see that the triviality of QED is in- 
timately connected with chiral symmetry breaking. Any 
attempt to remove the cut-off is always thwarted by the 
dynamically generated fermion mass. In particular this 
means that spinor QED does not exist as an interacting 
theory, similar to what Coleman and Weinberg |l(| found 
for scalar QED. 

We have also seen that chiral symmetry breaking al- 
lows QED to escape the Landau pole problem. While 
the bare parameters of the theory can take any values, 
the renormalized parameters are restricted. The Landau 
pole and ghost problem only occur deep in the inaccessi- 
ble e 2 R , am R region. Chiral symmetry breaking is always 
strong enough to push the Landau pole above the cut-off. 
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